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Abstract 

I summarize X-ray diagnostic studies of cosmic star 
formation in terms of evolutionary schemes for X-ray bi- 
nary evolution in normal galaxies with evolving star for- 
mation. Deep X-ray imaging studies by Chandra and XMM- 
Newton are beginning to constrain both the X-ray lumi- 
nosity evolution of galaxies and the logA^-log^ diagnos- 
tics of the X-ray background: I discuss these in the above 
context, summarizing current understanding and future 
prospects. 
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1. Introduction 

This is a brief account of the current status and future 
potentials of the X-ray diagnostics of the history of cos- 
mic star- formation rate (SFR). Global SFR has under- 
gone strong cosmological evolution: it was ~ 10 times its 
present value at z « 1, had a peak value ~ 10-100 times 
the present one in the redshift range z ~ 1.5-3.5, and 
declined again at high z (Madau, Pozzetti & Dickinson 
1998, henceforth M98; Blain, Smail, Ivison & Kneib 1999, 
henceforth B99a; Blain et al. 1999, henceforth B99b, and 
references therein). Details of the SFR at high redshifts 
are still somewhat uncertain, because much of the star 
formation at 2 <5 z ^ 5 may be dust-obscured and so 
missed by optical surveys, but detected readily through 
the copious submillimeter emission from the dust heated 
by star formation. 

The X-ray emission of a normal galaxy (i.e., one with- 
out an active nucleus) is believed to be dominated by the 
integrated emission of the galaxy's X-ray binary popula- 
tion: this statement may be somewhat dependent on the 
X-ray energy band (see Sec. ||), but current understanding 
does suggest that, in the canonical 2-10 keV X-ray band, 
the statement is valid for most normal galaxies. I sum- 
marize in this paper recent studies made in collaboration 
with N. White, A. Ptak, and R. Griffiths (White & Ghosh 
1998, henceforth WG98; Ghosh & White 2001, henceforth 
GWOl; Ptak et al. 2001, henceforth PtakOl) on the basic 
imprints of an evolving SFR on the evolution of X-ray bi- 
nary populations of galaxies, on the general consequences 



of these studies for deep X-ray imaging of galaxy fields 
by Chandra and XMM-Newton, and on the first, specific 
results that have emerged so far on the X-ray luminosity 
evolution in the Hubble Deep Field (HDF), and on the 
logiV-logS" diagnostics of the X-ray background. First re- 
sults of Brandt et al. (2001, henceforth BranOl) from the 
~ 0.5 Ms Chandra exposure of HDF North (HDF-N) sug- 
gest an evolution of the X-ray luminosities, Lx, of bright 
spirals from the Local Universe to z « 0.5, which I com- 
pare with the GWOl predictions from current SFR mod- 
els: I also discuss the roles of global and individual SFR 
histories in this context. Fluctuation analyses of the ~ 1 
Ms Chandra exposure of (HDF-N) suggest (Miyaji & Grif- 
fiths 2002, henceforth MG02) that the log TV-log S* plot in 
the soft X-ray band continues to rise at low {S ~ 10^^^ 
- 10~^^ erg cm~^ s~^) fluxes, indicating that the X-ray 
background at these fluxes is possibly dominated by a new 
population of faint X-ray sources rather than the canonical 
integrated AGN population (Gilli et al. 2001, henceforth 
GSH), whose contribution shows a cosmological flattening 
at these fluxes, and so lies much below the observational 
values suggested by MG02: in view of the PtakOl predic- 
tions which I discuss, it is plausible that this additional 
population is, in fact, that of normal galaxies showing the 
signature of their SFR histories through X-ray emission. 

2. SFR PROFILES AND X-RAY LUMINOSITY EVOLUTION 

In the approach of WG98 and GWOl, the total X-ray out- 
put of a normal galaxy is modeled as the sum of those of 
its high-mass X-ray binaries (HMXB) and low-mass X- 
ray binaries (LMXB), the evolution of each species "j" 
being described by a timescale . The effects of the depen- 
dence of Ti on the binary period and other parameters are 
studied by running the evolutionary scheme over ranges 
of likely values of given in the literature. The evolu- 
tion of the HMXB population in response to an evolving 
star-formation rate SFR(t) is given by 
9nHMXB(i) cT7TD/^^ "-hmxb(0 



ahSFRit) 



(1) 



dt thmxb 
where uhmxe is the number density of HMXBs in the 
galaxy, and thmxb is the HMXB evolution timescale. ah 
is the rate of formation of HMXBs per unit SFR, given 
approximately by an = 5/binary/piim/sN' where /binary 
is the fraction of all stars in binaries, f!}^:^ is that frac- 

' *' prim 

tion of primordial binaries which has the correct range of 
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stellar masses and orbital periods for producing HMXBs 
(van den Heuvel 1992, henceforth vdH92), and /gj^ w 1 is 
that fraction of massive binaries which survives the first 
supernova. In these calculations, a representative value 
THMXB 5 X 10^ yr is adopted according to current evo- 
lutionary models. Note that thmxb includes both (a) the 
time taken (~ 4 — 6 x 10^ yr) by the massive companion of 
the neutron star to evolve from the instant of the neutron- 
star-producing supernova to the instant when the "stan- 
dard" HMXB phase begins, and, (b) the (much shorter) 
duration (~ 2.5 x lO'^ yr) of this HMXB phase (vdH92 
and references therein). 

Of the two mechanisms of LMXB production generally 
envisaged, wz2;.,(a) production in cores of globular clusters 
due to tidal capture, and, (b) general production by evo- 
lution of primordial binaries, I describe here only the lat- 
ter one (which must be the dominant mechanism at least 
for spiral galaxies, since globular-cluster LMXB popula- 
tions of such galaxies can account only for relatively small 
fractions of their total X-ray luminosities), deferring the 
former to Sec. ||. LMXB evolution from primordial bina- 
ries has two stages (WG98) after the supernova produces 
a post-supernova binary (PSNB) containing the neutron 
star. First, the PSNB evolves on a timescale tpsnb due to 
nuclear evolution of the neutron star's low-mass compan- 
ion and/or decay of binary orbit due to gravitational ra- 
diation and magnetic braking, until the companion comes 
into Roche lobe contact and the LMXB turns on. Subse- 
quently, the LMXB evolves on a timescale tlmxb- Since 
TTSNB and tlmxb are comparable in general, the two stages 
are described separately (WG98) by: 

dnpsNBit) cT7T3/^\ "-PSNB(i) 



dt 

driLMXBit) 



aiSFR{t) - 
"psnb(^) 



TpSNB 

"-lmxb(0 



(2) 
(3) 



dt TpsNB Tlmxb 

Here, npgNB and tt-lmxb are the respective number den- 
sities of PSNB and LMXB in the galaxy, and ai is the 
rate of formation of LMXB per unit SFR, given approx- 
imately by ai = ^/binary/prim/sN' ^^c individual factors 
having meanings closely analogous to those for HMXBs 
(see GWOl). 

Evolution is displayed in terms of the redshift z, which 
is related to the cosmic time t by ig = 13(z -I- 1)"^/^, 
where tg is t in units of 10^ yr, and a value of Hq = 50 
km s^^ Mpc^^ has been used[|. I consider the suite of 
current SFR models detailed in Table 1 to cover a plausible 
range, using the parameterization of B99a,b. Models of the 
"peak" class have the form: 

-1 



(4) 



SFRpeak(^) = 2 1 + exp (l + 2)P+2l=r 

\ Zmax / 

^ For ease of comparison with WG98, M98, and GWOl, I use 
here a Friedman cosmology with go = 1/2. Other values of the 
Hubble constant lead to a straightforward scaling: for Ho = 70 
km s~^ Mpc~^, for example, tg ~ I0{z + 1)"^''^, so that the 
results remain unchanged if all timescales are shortened by a 
factor of 1.3. 



while those of the "anvil" class have the form: 

(1 ~t- ^)^, ^ ^ Zjjiax^ 

(1 ^" ^max)^ ^ ^ ^ Zmax- 



SFRa 



(5) 



Table 1. Star Formation Rate (SFR) Profiles 



Model 


Zmax 


P 


Comments 


Peak-M 


0.39 


4.6 


Madau profile 


Hierarchical 


0.73 


4.8 


Hierarchical Clustering 


Anvil-10 


1.49 


3.8 


Monolithic Models 


Peak-G 


0.63 


3.9 


Peak part of composite 








"Gaussian" Model 


Gaussian 


N/A 


N/A 


Gaussian starburst 








added at high z 



These functional forms are convenient since they have 
a convenient low-z limit, SFR(z) oc (l-l-z)^, where all SFR 
profiles must agree with the optical/UV data (M98), and 
since the model parameters can be manipulated to mimic 
a wide range of star- formation histories (B99b). Peak-class 
profiles are useful for describing (a) SFRs determined from 
optical/UV observations, i.e., Madau-type (M98) profiles, 
called "Peak-M" in Table 1, and, (b) more general SFRs 
with enhanced star formation at high z, an example of 
which is the "hierarchical" model of B99b, wherein the 
submillimeter emission is associated with galaxy mergers 
in an hierarchical clustering model. Anvil-class profiles are 
useful for describing the results of "monolithic" models. 
The "Gaussian" model (B99a,b)is an attempt at giving a 
good account of the SFR at both low and high z by mak- 
ing a composite of the Peak-G model (see Table 1) and a 
Gaussian starburst at a high redshift Zp, i.e., a component 

[t{z)-t{zp)f 



SFRGauss(z) = 6exp 



2cr2 



(6) 



Based on the IRAS luminosity function, this componenet 
is devised to account for the high-z data, particularly the 
submillimeter observations (B99a). For its parameters (see 
Table 1), I have used the revised values given in B99b (see 
GWOl). In all models described here, no galaxies exist 
for sufficiently large redshifts, z > 10. Figure 1 shows the 
prompt evolution of HMXBs and the slow evolution of 
LMXBs, and the evolution of the total X-ray binary pop- 
ulation, where the two components have been so weighted 
as to represent the total X-ray emission from the galaxy. 
The HMXB profile closely follows the SFR profile because 
Thmxb is small compared to the SFR evolution timescale. 
By contrast, the LMXB profile has a significant lag behind 
the SFR profile because tpsnb and tlmxb are comparable 
to SFR evolution timescale: the LMXB profile generally 
peaks at redshifts ^1-3 later than the HMXB profile — 
a charcteristic signature of SFR evolution ( WG98) . Effects 
of both (a) varying the evolutionary timescales for fixed 
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Table 2. Evolution of X-ray Luminosity Lx 
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Figure 1. Evolution of HMXB population (dotted line), 
LMXB population (dash-dotted line), and the total X-ray 
luminosity Lx (solid line) of a galaxy with various SFR 
profiles (dashed line), from GWOl. The effects of SFR 
variation are shown by keeping the evolutionary timescales 
fixed at tpsnb = L9 Gyr and tlmxb — 1-0 Gyr for all 
cases, and choosing various SFR profiles from Table L 
Each panel is labeled by the name of its SFR profile. 



SFR profiles, and, (b) varying the SFR profile for fixed 
evolutionary timescales have been studied: see GWOl for 
details. I display the latter variation in Figure 1 to em- 
phasize that, since, for sufficiently slow LMXB evolution, 
the galaxy's X-ray emission is dominated by LMXBs at 
low redshifts (0 ^ z ^ 1), and by HMXBs at high red- 
shifts, the total Lx-profile is strongly influenced at high 
redshifts by the SFR profile. Thus, determination of the 
Ljf -profile even upto moderate redshifts may put interest- 
ing constraints on the SFR, making this an independent 
X-ray probe of cosmic star-formation history. 

From their stacking analysis (see BranOl and refer- 
ences therein for an exposition of the technique), BranOl 
estimate that the average X-ray luminosity of the bright 
spiral galaxies at an average redshift z w 0.5 used in their 
study is about a factor of 3 higher than that in the local 
Universe. This observed evolution, Lx(0.5)/Lx(0.0) ^ 3, 
can be compared with the theoretical results in Table 2. 
The degree of evolution from z = to z = 0.5-1.0 in- 
creases from Madau-type profiles to those with additional 
star formation at high redshifts, the numbers for the Peak- 
M profile being in best agreement with BranOl. 



3. Global and individual SFRs 

A new development in SFR research in the last three years 
has been the study of star-formation histories of individ- 
ual galaxies and various galaxy-types. SFR profiles of in- 
dividual galaxies, ranging from those in the Local Group 
to those in the HDF at redshifts 0.4 ^ z ^ 1, have been 



Model 


TpSNB 


Tlmxb 


Lx(0.5) 
Lx(O.O) 


Lx(l.O) 
ix(O.O) 


Peak-M 


1.9 


0.1 


3.9 


5.4 


Peak-M 


0.9 


0.5 


4.6 


6.8 


Peak-M 


1.9 


1.0 


3.4 


4.1 


Hierarchical 


1.9 


1.0 


6.2 


11.3 


AnviI-10 


1.9 


1.0 


5.4 


10.1 


Gaussian 


1.9 


1.0 


7.5 


16.0 



inferred, using a variety of techniques. For various galaxy- 
types, models of spectrophotometric evolution, which use 
the synthesis code Pegase and are constrained by deep 
galaxy counts, have been developed (Rocca-Volmerange 
and Fioc 2000,henceforth RFOO), leading to a model SFR 
profile for each type. In the light of these developments, 
let me now suggest what may be the true significance of 
the BranOl results discussed above. 

BranOl used bright spirals for their stacking analysis. 
RFOO have shown that the model SFR profile for such (Sa- 
Sbc) spirals rises roughly in a Madau fashion from z = to 
z « 1 (which these authors ascribe to a bias in the original 
sample used to construct the Madau profile towards bright 
spirals), and thereafter flattens to a roughly constant value 
~ 12 times that at z = 0, falling again at z > 7. In the 
range < z ^ 7, this proflle can be roughly represented by 
an anvil- type profile (see Sec. |^), with the parameter Zmax 
as given in Table 1, and the parameter p k, 2.7. For such 
a profile with the timescales tpsnb = 1-9 Gyr, tlmxb = 
1.0 Gyr, as in Figure 1, the GWOl evolutionary scheme 
gives Ljf (0.5)/Ljf (0.0) = 3.3, and ix(l-0)/Ljf (0.0) = 5.4, 
in good agreement with both the BranOl results and the 
Peak-M results given in Table 2. It is now easy to see 
to see why the Peak-M profile would appear to give a 
good account of the BranOl results. In effect, the BranOl 
analysis may be probing the SFR profile of only the bright 
spirals in HDF-N, and the fact that the Peak-M profile is 
consistent with the BranOl results does not imply that the 
global SFR necessarily follows the Peak-M profile. 

4. log TV-log 5* DIAGNOSTICS: X-RAY BACKGROUND 

Based on the results of Sec. ^ PtakOl calculated the X-ray 
flux distributions and source count (log A^-log S) plots ex- 
pected for HDF-N. Figure 2 shows the PtakOl plot in the 
soft (0.5-2.0 keV) X-ray band, which has proved to be a 
valuable diagnostic of current population synthesis models 
of the X-ray background, as I now summarize. Hasinger 
(2002, henceforth II02) reminds us that the cosmic X-ray 
background has, in this Chandra and XMM-Newton era 
of deep X-ray surveys, been largely resolved into contri- 
butions from individual sources, the resolved fraction be- 
ing ^ 90% in the soft (0.5-2.0 keV) band, and similar 
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Figure 2. log -log 5 plot in the soft (0.5 - 2.0 keV) band 
for HDF-N, from PtakOl. The diamonds correspond to the 
Gaussian SFR profile described in Sec. and the crosses 
to the Peak-M profile. Note that an interpolation through 
the former points is represented by a dashed line in Fig. 2 
of MG02 and Fig. 2 (bottom) of H02, and that through 
the latter points by a dotted line. The solid line here is 
the double power law fit of Tozzi et al. to the Chandra 
observations of HDF South. 



in the harder (2-10 keV) band. The long-standing belief 
that these sources are predominantly active galactic nu- 
clei (AGN), both unobscured and obscured (the so-called 
QS0-2s), was supported by the (now completed) optical 
identification programme which followed up the ROSAT 
deep survey, since it found the counterparts to be predom- 
inantly AGN. Ongoing optical identifications of the deep- 
est Chandra and XMM-Newton fields are still far from 
complete. AGN population-synthesis models of the X-ray 
bacground are currently very useful and popular: these 
have been developed to a degree of sophistication and de- 
tail (see GSH, which has references to earlier models) suf- 
ficient for extracting information about AGN population 
properties. The recent, ultradeep (~ 1 Ms) observations of 
both HDF-N and the Chandra Deep Field South (CDFS) 
have led to log TV-log S" plots in the soft (0.5-2.0 keV) X- 
ray band which go down to fluxes ^ 5 x 10^^^ erg cm~^ 
s""'^: these are fitted well by the GSH models, which show 
a clear cosmological flattening at fluxes below the above 
limit (H02; MG02). 

Fluctuation analysis is a powerful tool for constraining 
the source counts below source detection limit (see MG02 
and references therein for an exposition of the method), 
which has been successfully tested on data from previ- 
ous X-ray missions. Its recent application by MG02 to the 
1 Ms observation of HDF-N has yielded the remarkable 
result that the constraints so obtained on the soft-band 
logA^-logS plot suggest that the extension of the plot 
down to fluxes 5* ~ 7 x 10~^* erg cm~^ s~^ continues to 
rise as at higher fluxes, showing no signs of the cosmolog- 
ical flattening characteristic of the GSH models (MG02; 



H02). The most obvious interpretaion is that, while the 
AGN contribution, as modelled by GSH, begins to sat- 
urate at these fluxes, a new population of faint sources 
begins to dominate. The fact that the the extension of the 
logiV-logS" plot, as inferred from the fluctuation- analysis 
constraints of MG02, agree well with that shown in the 
above PtakOl plot (since the figure showing this appears 
twice in these proceedings. Fig. 2 in MG02 and Fig. 2 (bot- 
tom) in H02, I do not repeat it here), particularly for the 
Gaussian SFR profile, therfore opens the exciting possibil- 
ity that first signatures of cosmic star formation in the soft 
X-ray band logiV-logS' plots are revealing themselves. 

5. Future prospects 

Even in this Chandra and XMM-Newton era, truly new 
results on Lx-evolution and SFR signature have been pos- 
sible so far only by going below the source detection limit 
with special techniques like stacking and fluctuation anal- 
ysis. These are suggestive indications, which must be con- 
firmed with source detection at lower fiuxes, first with 
longer exposures with Chandra and XMM-Newton, and 
then with the next generation of satellites like Constella- 
tion-X and XEUS. On the theoretical side, the evolution- 
ary scheme must be generalized to include several addi- 
tional effects, e.g., (a) in the soft X-ray band, the output 
of a normal galaxy may have very significant contributions 
from supernova remnants (G. Hasinger, personal commu- 
nication), and, (b) tidal capture creation of LMXBs in 
globular clusters may be the dominant production mech- 
anism in certain galaxy-types. Inclusion of these effects 
presents no difficulties of principle, and is now under way: 
the results will be described elsewhere. 
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